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Abstract. TeV gamma-rays were recently detected from the nearby normal spiral galaxy NGC 253 (Itoh et al., 
2002). Observations to detect the Cherenkov light images initiated by gamma-rays from the direction of NGC 253 
. ^ were carried out in 2000 and 2001 over a total period of ~150 hours. The orientation of images in gamma-ray-like 

^ . events is not consistent with emission from a point source, and the emission region corresponds to a size greater 

than 10 kpc in radius. Here, detailed descriptions of the analysis procedures and techniques are given. 

Key words. Gamma rays: observation - Galaxies; starburst - Galaxies: individual: NGC 253 - Galaxies: halos - 
cosmic rays 
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1. Introduction 

NGC 253 is a very nearby {d = 2.5 Mpc) (de Vaucouleurs, 
I1978|l . normal spiral, starburst galaxy. Starburst galaxies 
are generally expected to have cosmic-ray energy densi- 
ties about hundred times larger than that of our Galaxy 
(Voelk et al.. 11989(1 due to the high rates of massive star 
formation and supernova explosions in their nuclear re- 
gions. The star-formation rates can be estimated from the 
far- infrared (FIR) luminosities, and the supernova rates 
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can be also inferred based on the assumption of an ini- 
tial mass function. Since the supernova rate of NGC 253 
is estimated to be about 0.05 - 0.2 yr~^ (Mattila and 
Meikle, IMTTl Antonucci and Ulvestad, HMSl van Buren 
and Greenhouse. II99H1 . a high cosmic-ray production rate 
is expected in this galaxy. 

We recently reported on the detection of TeV gamma- 
rays from NGC 253 (Itoh et al. 2002). Previous to this, 
the only evidence for higher energy particles in a galaxy 
other than our own is for the Large Magellanic Cloud 
(Sreekumar et al., I1992|l . Voelk et al. H1996fl estimated 
the gamma-ray fluxes (via neutral pion decay) from the 
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nucleus of nearby starburst galaxies. These values, how- 
ever, were under the sensitivity of the EGRET detector 
on the Compton Gamma-Ray Observatory (CGRO) and, 
indeed, EGRET observations resulted in very stringent 
upper limits for the GeV emission from NGC 253 (Blom 
et al. lTMl Sreekumar et aL. ITUMjl . 

On the other hand, NGC 253 has an extended 
synchrotron-emitting halo of relativistic electrons (Carilli 
et al., I1992|l . The halo extends to a large-scale height, 
where inverse Compton scattering (ICS) may be a more 
important process for gamma-ray production than pion 
decay and bremsstrahlung. The seed photons for ICS are 
expected to be mainly FIR photons up to a few kpc from 
the nucleus, and cosmic microwave background radiation 
at larger distances. 

The OSSE instrument onboard the CGRO detected 
sub-MeV gamma-rays from NGC 253 (Bhattacharya et 
al., '1994!). This emission is consistent with a model 
for ICS of the FIR photons around the nucleus of the 
galaxy by synchrotron-emitting electrons (Goldshmidt 
and Rephaeli. I1995|l . although it is difficult to study the 
spatial distribution of the emission due to the limited an- 
gular resolution of OSSE. 

We observed NGC 253 with the CANGAROO-II tele- 
scope in 2000 and 2001, and detected TeV gamma-ray 
emission at high statistical significance (Itoh et al. 2002). 
This detection of TeV gamma-rays from a normal spiral 
galaxy like our own has profound implications for the ori- 
gin and distribution of TeV cosmic-rays in our Galaxy. 
In this paper we describe in detail the observations and 
analysis of the TeV gamma-rays from NGC 253. 

2. Observations 

2.1. The CANGAROO-II telescope 

The CANGAROO (Collaboration of Austraha and 
Nippon (Japan) for a GAmma Ray Observatory in 
the Outback) air Cherenkov telescope is located near 
Woomera, South Australia (136°46'E, 31°06'S, 220 m 
a.s.l.). The telescope consists of a 10 m refiector and a 552 
pixel camera. It detects images of cascade showers result- 
ing from sub- TeV gamma-rays (and background cosmic 
rays) interacting with the Earth's upper atmosphere. 

The CANGAROO-II project is exploring the south- 
ern sky at gamma-ray energies of 0.3^100 TeV. Its pre- 
decessor, CANGAROO-I, used a 3.8 m telescope (Hara et 
al., HM^ . and detected TeV gamma-ray emission from 
such objects as pulsar nebulae (PSR 1706-44 (Kifune 
et al. I1995|l . the Crab (Tanimori et al.. I1998a|l '). super- 
nova remnants (SNR) (SN1006 (Tanimori et al.. I1998b|l . 
and RXJ1713. 7-3946 (Muraishi et al.. lMin| ). The 10m 
telescope of CANGAROO-II has been in operation since 
April, 2000, a nd has detected SNR RXJ1713. 7-3946 
(Enomoto et al. l2002b|l and the active galactic nuclei Mrk 
421 (Okumura et al. . 1200211 . The telescope has a parabolic 
optical reflector consisting of 114 composite spherical mir- 
rors (80cm in diameter), made of carbon fiber reinforced 



plastic (CFRP) (Kawachi et al.. l2001|l . The principal pa- 
rameters of the telescope are listed in Table 



Parameters 



Values 



Location 136 °E ,31°S 

Height above sea level 220 m 

Total diameter 10 m 

Focal length 8 m 

Number of segmented mirrors 114 

Mirror diameter 80 cm 

Mirror segment shape Spherical 

Mirror alignment Parabolic 

Mirror curvature 16.4 m 

Mirror material Plastic (CFRP) 

Table 1. Principal parameters of the CANGAROO-II 
telescope. 



The camera contains 552 pixels, each of which sub- 
tends an areaof 0.115°x0.115°. Each pixel is a 1/2" photo- 
muhipher tube (PMT) (Hamamatsu Photonics R4124UV) 
with an air light guide. The output signal is amplified by 
a high-speed IC (Lecroy TRA402S) and split three ways 
for the ADC (analogue to digital converter), TDC (time 
to digital converter), and the scalers. The scaler is a spe- 
cial circuit which records the number of hits greater than 
the threshold (> 2.5 photoelectrons) of individual PMTs 
within 700/isec (Kubo et al. l2001|l . The scalers were trig- 
gered by a clock (1 Hz), and these data were recorded 
every second. 

2.2. Pointing direction and observation 

The telescope was pointed at the center of 
NGC 253, the J2000 coordinates of which are (RA, 
Dec) = (11.888°, -25.288°). NGC 253 was observed from 
October 3 to November 18, 2000 and from September 20 
to November 15, 2001, with the CANGAROO-II tele- 
scope. The observations were carried out on clear nights 
during moon- less periods. Periods of 1.5 hours after sun- 
set and 1.5 hours before sunrise were avoided. Each night 
was divided into two or three periods, i.e., ON-OFF, 
OFF-ON-OFF, or OFF-ON observations. ON-source 
observations were timed to contain the meridian passage 
of NGC 253, which culminates at a zenith angle of 6°. 
The observation times are summarized in Table |21 In 



Observation Date 


Ton (min) 


To// (min) 


03-Oct. - 18-Nov. 2000 


2297 


2245 


20-Sep. - 15-Nov. 2001 


2567 


2401 


Total 


4864 


4646 



Table 2. Summary of the observation periods. 



total, ~4800 minutes of ON-source observations and a 
similar amount of OFF-source observations were carried 
out. 
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2.3. Hardware Trigger 

The pixel arrangement of the CANGAROO-II camera is 
shown in Fig.^ The trigger-region is defined by the inner 

CAHGAROO-II CAMERA 36 Boxes 




Fig. 1. Pixel arrangement of the CANGAROO-II camera. 
The thick Box (1.84°x 1.84°) is the trigger region. The 
camera consists of 36 boxes. Each box contains 16 PMTs, 
each 1/2" in diameter. In total, 552 PMTs are installed. 
Each pixel subtends 0.115°x0.115°, defined by the light 
guide. 



photoelectron, the cosmic-ray rate roughly agreed. Using 
a Monte-Carlo simulation (described later) of cosmic-ray 
protons, we further studied this ADC conversion factor 
fHara. l?()()2|l . We analyzed the relation between the total 
ADC counts and the total number of pixel hits. From this 
correlation we determined this factor to be 92 l^'^ [ADC 
ch/p.e.]. This agreed with the results of a study of the 
Night Sky Background (NSB) rate. 



3.3. Pixel Selection 

Occasionally, individual pixels display anomalously high 
count rates. The trigger rate of each pixel is monitored 
by a scaler every second. This information enabled us 
to remove these "hot" pixels from further analysis. The 
scaler distributions obtained from 2000- and 2001-data are 
shown in Figs. |21a) and b), respectively. In 2000, due to 



1.84°x 1.84° square, which contains 256 PMTs in 16 boxes. 
The event trigger requires: 

1. More than three pixels to be hit inside the trigger re- 
gion. The threshold for each pixel was set at approxi- 
mately 2.5 photoelectrons (p.e.). 

2. More than one box with a charge-sum exceeding 
-10 (p.e.). 

3. Data analysis 

3.1. Calibration 

The data were calibrated using a LED (Light Emitting 
Diode) light source located at the center of the 10 m mir- 
ror, ^--^8 m from the camera (Kabuki et al., .2002). A 
quantum-weU type blue LED (NSPB510S, A ~ 470 nm, 
Nichia Corporation, Japan) was used, and illuminated 
with an input pulse of ^ 20 nsec width. A light diffuser 
was placed in front of the LED in order to obtain a uni- 
form yield on the focal plane. The main purpose of this 
calibration was field flattening. The relative gain of each 
pixel was adjusted according to the mean pulse height of 
all pixels. The second purpose was to adjust the timing of 
each pixel with respect to the mean timing for all pixels. 
Time-walk corrections (adjusting the earlier triggering of 
larger pulses that arises from a fixed trigger threshold) 
were also carried out, based on the data. This calibration 
was done run by run. 

3.2. ADC conversion factor 

In order to compare the simulated and observed spectra, 
the energy scale must be calibrated, i.e. a conversion fac- 
tor from the ADC value to the absolute energy is required. 
First, we checked the cosmic ray event rate. Under the as- 
sumption that ~ 100 ADC counts corresponded to a single 




10 20 10 20 

Scaler Scaler 



Fig. 2. Scaler distribution: a) 2000-data and b) 2001-data, 
where "N" denotes the number of hits per scaler count 
(horizontal axis). The rate at which each pixel exceeds 
the ~ 2.5 p.e. threshold was monitored for 700 /iS each 
second. 



the influence of artificial lighting from the detention cen- 
tre several km away, the hit rate was significantly higher 
than that of 2001, and so a slightly higher cut value was 
adopted for the 2000 data. There were no stars brighter 
than a magnitude of 5.6 in the field of view (FOV) of the 
camera during these observations. However, the effects of 
fainter stars passing through the FOV of a pixel were ex- 
pected to be removed by this PMT rate cut. This was 
confirmed using data from other observations which had 
brighter stars in the FOV of the camera. 
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3.4. Clustering 



The purposes of the pre-selection were to remove noisy 
pixels affected by the NSB and any period affected by 
cloudy conditions from the observation data. Here, we 
used "tna" logic (threshold n-adjacent, where n is the 
number of adjacent PMTs required to have triggered) . The 
threshold was fixed at around 300 ADC count (approxi- 
mately 3.3 p.e.). The distribution of ADC is shown was 
in Fig. |31 The hardware threshold was located at ~200 
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Fig. 3. ADC distribution for pixels, after pedestal levels 
had been subtracted, where "N" denotes the number of 
hits per 10- ADC count. The small peak in negative region 
is due to the electronics undershoot. 

ADC counts. Therefore, the cut at 300 ADC counts (~3.3 
p.e.) is reasonable. This is the simplest and most powerful 
method to reject pixels affected by NSB. 

After this selection, those clusters with more than n 
adjacent hits were selected. As n increased, the TDC dis- 
tribution became cleaner, as shown in Fig. ^ The mean 
event timing was located at around 300 TDC counts 
(1 count = 1 nsec). Those events uniformly distributed 
between 200 and 400 nsec are considered to be due to 
NSB photons. From Fig. 01 we selected a cut of n = 4. We 
also cut pixels with |TDC-300| > 40 nsec. 

After this cut, the event rate which satisfies t4a-cluster 
and does not satisfy t5a is shown in Fig. |31c). Although 
the raw trigger rates were not stable, due to changes in the 
background light level as the telescope pointing changed, 
the shower rate became stable when t4a-selection was ap- 
phed, as can be seen in Fig. Elc). The shower rate after 




100 200 300 400 500 

TDC distribution (nsec) 



Fig. 4. TDC distributions: a) for all pixels, b) for pixels 
which satisfy t3a-clustering and do not satisfy t4a, and c) 
for pixels which satisfy t4a-clustering and do not satisfy 
t5a. "N" denotes the number of hits per TDC count, which 
were normalized to 1 nsec. The horizontal axis is the TDC 
count. 

the t4a-clustering is shown in Fig. Eld). From then on, the 
events with t4a-clustering were selected. 

Using these shower rate plots, we were also able to re- 
move any cloudy period during an observation. Examples 
of these plots for good and bad conditions are shown in 
Figs. El a) and c), respectively. The cloudy periods de- 
tected in this manner perfectly matched the observing 
conditions described in the experimental log. The cut line 
(the dashed line in Fig. jSJ was set at 2.0 Hz, correspond- 
ing to 600 events per five minutes. The shower rate of 
the data passing these cuts was very stable over all ob- 
servations, in all seasons, except for the expected zenith 
angle dependence. The stability of the shower rate in 2001 
is shown in Fig. jJl This stability implies that the degra- 
dation of the mirror reflectivity was small. Over a short 
period in 2000, the camera system was exposed to strong 
background light from the township of Woomera, and the 
shower rate was observed to change by ~20%. This cor- 
responds to a change in the energy scale of 10% for the 
i?^^'^ cosmic-ray spectrum. A systematic error of 10% is 
included in the energy determination as a result. Also, 
we rejected any events with zenith angles greater than 20 
degrees, as shown in Figs. Ob) and d). 

These observations were carried out in the south- 
ern hemisphere during spring. Near sunrise, the humid- 
ity increased, and dew started to form on the surfaces of 
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Fig. 5. Shower rates: a) hardware trigger rate; b) shower 
rate for events satisfying t3a-clustering, but not satisfying 
t4a; c) shower rate for events satisfying t4a-clustering, but 
not satisfying t5a; and d) shower rate for t4a-clustering. 
The vertical axis " N" denotes the number of events per 5 
min. The horizontal axis is the time in minutes from the 
start of the observation. 



the mirrors. These effects could also be detected by this 
shower-rate study, and we eliminated these periods from 
any further analysis. 

After these pre-selections, the data remaining for anal- 
ysis were accumulated, as summarized in Table |3| 



Observation period 


Ton (min) 


To// (min) 


Ton/To// 


3 Oct. - 18 Nov. 


1301 


969 


1.34 


2000 








20 Sep. - 15 Nov. 


1658 


1448 


1.15 


2001 








Total 


2959 


2417 


1.22 



Table 3. Summary of data remaining after pre-selection 
cuts. 
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Fig. 6. Shower rate under a) good conditions and c) poor 
conditions. The horizontal axis is the time (minutes) from 
the start of an observation. The vertical axis ("iV") is 
the number of events per 5 minutes. The histogram is the 
shower rate, and the dashed line is the cut position, i.e., 
600 events per 5 minutes were required for the data to be 
accepted. Plots b) and d) are the elevation angle distribu- 
tions versus time. Data was accepted above an elevation 
angle of 70°. 

3.5. Image Analysis 
3.5.1. Monte-Carlo Setups 

Simulations of electromagnetic and hadronic showers in 
the atmosphere were carried out using a Monte-Carlo sim- 
ulation code based on GEANT3.21 HGEANT|I . In this 
code, the atmosphere was divided into 80 layers of equal 
thickness (~ 12.9g/cm^) (Enomoto ct al., .2002afl . Each 
layer corresponds to less than a half radiation length. 
The dependence of the results on the number of layers 
was checked by increasing the number of layers, and was 
confirmed to be less than a 10% effect. The lower energy 
threshold for particle transport was set at 20 MeV, which 
is less than the Cherenkov threshold of electrons at ground 
level. Most Cherenkov photons are emitted higher in the 
atmosphere, at lower pressure and a higher Cherenkov 
threshold. The geomagnetic field at the Woomera site 
was included in the simulations (a vertical component of 
0.520 G and a horizontal component of 0.253 G directed 
6.8° east of south). 

In order to save CPU time, Cherenkov photons were 
tracked in the simulations only when they were initially 
directed to the mirror area. The average measured re- 
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Fig. 7. Shower rates versus the cosine of the zenith angle. 

flectivity of 80% at 400 nm and its wavelength depen- 
dence (Kawachi et al., 12001(1 and the measured PMT 
quantum efficiency were multiplied using the Frank- Tamm 
equation to derive the total amount of light and its 
wavelength dependence. A Rayleigh-scattering length of 
2970(A/400nm)(g/cm^) (Baum and Dunkelman, [TU33)l 
was used in transport to the ground. No Mie scattering 
was included in this study. The contribution of Mie scat- 
tering is thought to be greatest at the 10-20% level; we 
therefore consider this study to have uncertainties of at 
least this level. When Rayleigh scattering occurred, we 
treated it as absorption. 

Finally, the simulated electronic noise was added and 
the timing responses were smeared using a Gaussian of 
4 ns (Icr). We also added NSB photons, conservatively se- 
lecting to double Jelley's value of 2.55 x 10~*erg/cm^/s/sr 
(430-550 nm) (Jelley, I1958|l . Electronics saturation was 
also taken into account. The zenith angle distribution was 
obtained from an ON-source run, and is shown in Fig. |H1 
The above distribution was input to the event genera- 
tor of the Monte-Carlo simulation. We generated gamma- 
rays between 100 GeV and 10 TeV assuming a (Crab-like) 
power-law spectral index of —2.5. 

The energy spectrum for simulated events passing the 
pre-selection cuts are shown in Fig.|^ From this figure, we 
obtained the threshold energy for the gamma-ray detec- 
tion to be ~ 500 GeV for a iJ-^-S spectrum and 400 GeV 
for a E~^ '^ spectrum, after pre-selection cuts. These simu- 
lated gamma-ray events were used together with observed 
events from OFF-source runs to determine the cut values 
in order to optimize the gamma-ray signal. 



Fig. 8. Zenith-angle distribution after event selection. The 
vertical axis is the observation time in minutes. The mean 
is 11.1°. 

3.5.2. Data analysis 

We first calculated the standard image parameters: 
Distance, Width, and Length (Hillas, I1985|l . The distri- 
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Fig. 9. Distribution of energies for accepted events from 
the Monte-Carlo gamma-ray simulation; the solid line was 
obtained for a E~^-^ spectrum and the dashed line for a 
g-3.0 spectrum. 
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Fig. 10. Image parameter distributions. The blank his- 
tograms were obtained from OFF-source observations. 
The hatched areas are the distributions for gamma-rays 
from our Monte-Carlo simulations. 



to check the background shapes of image parameters, we 
carried out a simulation of cosmic-ray proton events, gen- 
erating protons (only) between 500 GeV and 10 TeV from 
a differential i?"^ '' spectrum. The mean elevation angle 
of OFF-source observations was assumed. The distribu- 
tions of the resulting image parameters were checked and 
found to be roughly consistent with those obtained by the 
OFF-source runs. 

We cut events with a Distance of less than 0.5° or 
greater than 1.2°. The Width and Length were used as 
a Likelihood ratio, which is described later. We defined 
Eratio as the sum of the ADC counts outside the main 
cluster in the image, divided by the sum of the ADC sum 
inside the main (maximum energy) cluster. Gamma-ray 
events are predicted by simulations to be typically a sin- 
gle cluster, and thus have low values of Eratio- The distri- 
bution of Eratio is shown in Fig. We rejected events 
with Eratio > 0.1. This cut helped to reduce the cosmic- 
ray backgrounds. Also, less-energetic events, with a total 
ADC count of less than 3000 (~33 p.e.), were rejected in 
order to improve the a (image orientation angle) resolu- 
tion. 

The acceptance of gamma-ray-like events was evalu- 
ated using the Likelihood-ratio (Enomoto et al., I2002al 
and l200l')l . Probability Density Functions (PDFs) were de- 
rived for both gamma-ray and cosmic ray initiated events. 
The PDFs for gamma-rays were obtained from simula- 
tions, while those for cosmic rays were obtained from OFF- 
source data. Histograms were made of Length and Width 
using both data sets; these distributions were then nor- 
malized to unity. The probability (L) for each assump- 
tion was thus obtained by multiplying ViyF^Width) by 
VIW [Length). In order to obtain a single parameter, and 
also to normalize it to unity, we used the Likelihood-ratio: 




Likelihood — ratio ~ 



L(gamma — ray) 



L{gamma — ray) -\- L{proton) 



(1) 



ratio 



Fig. 11. Eratio- Ratio between the ADC counts of the 
maximum energy cluster and that of the remaining clus- 
ters. "N" denotes the number of events. The blank his- 
togram was obtained from an OFF-source run and the 
hatched histogram from the gamma-ray simulations. 

We took care to take account of the energy depen- 
dences of these two image parameters. An example of the 
Width parameter is shown in Fig. ^] In order to cor- 
rect for the energy dependence, we made 2D-histograms, 
i.e., the shape parameter versus the logarithm of the to- 
tal ADC counts, and calculated PDFs for an energy in- 
dependent acceptance of signal events. Here, we did not 
use either the Distance or Asymmetry parameter, as these 
parameters are source-point dependent. If the gamma-ray 
emitting region is broader than that of a point source, 
these parameters will deviate from that of point source. 
In the Monte-Carlo simulations for the PDF determina- 
tion, we used the point-source assumption. 

We selected gamma-ray-like events using the 
Likelihood-ratio. The Likelihood-ratio distributions are 
shown in Fig. ^| The signal peaks at a Likelihood-ratio 
of 1 and the background at 0. 

We then investigated the figure of merit (FOM) using 
these data in order to maximize the statistical significance 
of the gamma-ray signals from NGC 253. At various cut 
locations, the signal of the Monte-Carlo simulation and 
the OFF background entry were obtained. The FOM was 
defined as the former value divided by the square root of 
the latter value. The FOM versus Likelihood-ratio cut val- 
ues are plotted in Fig.^^ The figure suggests that higher 
cut values lead to a higher statistical significance, albeit 
with a loss in the gamma-ray acceptance. As a compro- 



8 



Itoh, Enomoto, Yanagita, Yoshida et al.: TeV gamma-ray from NGC 253 




Fig. 12. Correlation between the Width and the logarithm 
of total ADC counts obtained from the gamma-ray simu- 
lations. 



mise between the acceptance and the FOM, we opted to 
adopt a value of 0.4 for the Likelihood-ratio cut, noting 
that there was only a small change in FOM between 0.2 
and 0.6. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Likelihood-ratio Likelifiood-ratio 



Fig. 13. Likelihood Distribution: a) for 2000 data and b) 
for 2001 data. "N" denotes the number of events. The 
blank histograms are OFF-source data and the hatched 
areas are from the gamma-ray simulations. 



3.6. Further Hot Pixel Rejection 

Pixels occasionally have anomalously high trigger rates, 
often due to enhanced starlight or man-made light in the 
FOV of the pixel, or to small discharges between the light- 
guide and the photo cathode (Kabuki et al. , '2002") , or to 
electrical noise in the associated circuitry. Although these 
are generally random, small pulse-height signals, the high 
pixel trigger rates can have the affect of increasing the 
camera trigger rate. When randomly triggered during a 
real event trigger, these "hot pixels" rarely form a clus- 
ter. Their effect was thought to be reduced after the pre- 
selection clustering cuts. However, it is possible that out- 
lying hot pixels surviving the pre-selection cuts deform 
the shapes of the shower images. Such effects could sig- 
nificantly smear the a distribution for gamma-ray events. 
In fact, the a distributions of the OFF-source runs were 
observed to be deformed from the Monte-Carlo prediction. 

In order to flatten these distributions, we removed hot 
pixels. First we looked at the hot pixel map for events 
passing selection cuts. This enabled the hottest pixels to 
be identified and removed. We then looked at the scaler 
counts for the remaining pixels. In the same way, some of 
the hottest pixels were removed. Finally, we tested pixels 
iteratively to find out whether they deformed the a dis- 
tributions of OFF-source run events for clusters having a 
center of gravity around the pixel being investigated. For 
the 2000 data, 12.3% of the pixels were removed by these 
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operations. The same procedure was carried out for the 
2001 data, and 9.7% of pixels were removed. 

After removing those hot pixels, we checked the Monte- 
Carlo simulation precisely, and verified that these proce- 
dures did not result in any deformation of the image pa- 
rameters, including a. 

Most of hot pixels were located around the edges of the 
camera. This occurred as PMTs with high trigger rates, 
which were identified early during camera testing, and 
were deliberately moved from inside the trigger region to 
the outer edge of the camera to minimize the effect on the 
hardware trigger. 

In order to check whether small pulse-height random- 
noise signals were removed by these operations, we loos- 
ened the clustering to t3a, which should be more sensitive 
to these backgrounds. Similar plots were obtained, which 
confirmed that they were still consistent with the Monte- 
Carlo predictions. With t4a clustering and these proce- 
dures, we concluded that the random noises were removed 
successfully by this operation. 

3.7. Results 

The resulting a distributions arc shown in Fig. 1151 The 



RX J1713.7— 3946, which was found to have an extended 
nature (Enomoto et al. I2002b)l . The distortion of the a 
spectrum in both ON- and OFF-source runs was observed 
in 2001 (Fig.H^b). The ON and OFF spectra, however, 
agreed well for a > 30°, even including the normalization 
factor, which is described in section ITTl These appeared 
in higher a regions, and were considered to be due to hot 
channels which remained even after the rejection proce- 
dure described in the previous section. Most of these hot 
channels were located outside of the trigger region. In or- 
der to keep the high efficiency of the analysis, the hot 
channels which had less affect on the deformation of the 
a spectrum were not removed. For a point source, simula- 
tions predict that an a cut of between 15 and 20° should 
optimize the signal. The simulations indicate that 73% 
of the excess from a point source should have a < 15°. 
We tested a cuts from 15 to 35° in 5° steps in order to 
maximize the excess for this source, and found that 30° 
was best. The statistical significance of the signal thus 
needs these 5 trials to be taken into account. The final 
significance remained greater than lOcr. The signal rates 
were (743±105) events/ 1301min = 0.57±0.08 for 2000 
and (933±106) event/ 1658min = 0.56±0.06 for 2001, con- 
sistent within the statistical errors. The average event rate 
overall was 0.56±0.05 /min. 




20 40 60 800 20 40 60 80 
a(degree) a{degree) 



Fig. 15. The distributions of a obtained for a) 2000 data 
and b) 2001 data. The points with error bars were ob- 
tained from ON-source data. The shaded histogram is 
from the OFF-source data. 

excesses at a < 30° are 742.5±104.6 events (7.1(7) for 
2000 and 933.1±106.2 events (8.8cr) for 2001, respectively. 
By combining the two years of data, the total excess 
was found to be 1651.9±149.2 (ll.lcr). Here, we used 
an Likelihood-ratio cut of 0.4. The a cut at 30° was 
larger than expected for a point source. This cut value 
was used in the previous detection of gamma-rays from 



3.8. Various Checks 

We investigated the effect of raising the Likelihood-ratio 
cut, 0.6, i.e., applying a tighter cut. The results for the 
combined (2000 and 2001) data set are shown in Fig. ^| 
b). The excess for the tight cut is 823.8±105.6 (7.8cr), 




20 40 60 800 20 40 60 80 20 40 60 80 



a(degree) a(degree) a(degree) 

Fig. 16. a distributions obtained with a) the loose cut 
(L > 0.4), b) a tighter cut {L > 0.6), and c) "standard" 
cuts. The points with error bars were obtained from ON- 
source data. The shaded histogram represents OFF-source 
data. 

somewhat lower, as expected due the reduced acceptance. 

In order to verify the Likelihood method, we checked 
the results of the "standard" analysis, using an accep- 
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tance cuts of 0.05< Width <0.15 and 0.05< Length <0.3. 
The a distributions are shown in Fig.^]c). The excess is 
1696.7±165.9 (10.2cr), with signals of 811.7±119.0 events 
(6.8cr) and 916.7±115.4 events (8. Oct) in 2000 and 2001, re- 
spectively. As expected, the standard analysis confirms the 
statistical significance of the detection, though at a lower 
level than the more powerful Likelihood-ratio method. 



Approximately 10 hours of ON- and OFF-data were used 
for analysis. The energy threshold was estimated from sim- 
ulations to be ~2 TeV. The excessed number of events 
was 405 ± 59 (6.8cr), as shown in Fig.UHla). The differ- 



3.9. Effective Area 

The effective areas for this analysis is shown in Fig. El We 
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Fig. 17. Effective area versus energy; a) Effective areas 
after the pre-selection (the black squares), those after 
the distance cut (the black triangles), those after the 
Likelihood-ratio cut (the blank circles), and those with 
a < 15° (the blank squares). The effective area for the 
Whipple telescope is shown by the line (after the distance 
cut), b) The effective areas multiplied by E~^-^ are shown 
in order to indicate the threshold of the CANGAROO 
telescope. 
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compared them with that of Whipple (Fig. 6 in Mohanty 
et al. 1998), which are the effective areas after clustering 
and Distance cut, i.e., before image parameter cut. Our 
effective area agreed with Whipple, even with the energy 
dependences. According to Fegan (|1996|l . the threshold of 
Whipple telescope was the same as ours. 

4. Discussion 

4.1. Crab Analysis 

We observed the Crab nebula in November and December, 
2000, in order to check our energy and flux determina- 
tion. The Crab nebula has a power-law spectrum over a 
wide energy range (Aharonian et al., '2UUU', Tanimori et 
al. . 11998611 . The elevation angle ranged from 34° to 37°. 



Fig. 18. Results of the Crab analysis. The differential 
flux obtained for the Crab nebula is shown in b), together 
with previous observations. The points with error bars 
were obtained by this experiment. The dotted line is the 
HEGRA resuh and the dashed hue is the CANGAROO- 
I result. The insert (a) is the a distribution. The points 
with error bars are obtained from ON-source data and the 
histogram is from OFF-source data. The lower plot (c) is 
a signiflcance map. The 65%-contour is drawn with the 
our estimated angular resolution (the arrows: ±1ct). 

ential flux, shown in Fig. 1181 -b'). is consistent with results 
from HEGRA (Aharonian et al.. . 2000,1 and CANGAROO- 
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I (Tanimori et al.. ll998a|l observations. We conclude that 
our energy and flux estimations are correct. In addition, 
we derived a cosmic-ray spectrum from background events 
and compared it with the known cosmic-ray flux. From 
these checks, the systematic uncertainty for the absolute 
flux estimation was found to be within 10%. 

The a plot for the Crab nebula, shown as an insert 
(Fig. ^la), is consistent with the point-source assump- 
tion. The OFF-a spectrum was again not flat. The average 
image positions were different from the zenith-angle obser- 
vation, i.e., they were centerized because the shower max 
position was higher in altitude. The reflective index of air 
was smaller there, which resulted in a smaller Cherenkov 
angle. In addition, there is a bright star close to the Crab 
position. In order to avoid a high trigger rate, we dis- 
placed (0.25 degrees) the telescope's tracking center away 
from the center position of the Crab nebula. The average 
hit region in the camera was different from the observa- 
tion of NGC 253. Thus, a different a deformation (due to 
the hot channels) occurred in this case. The 65% contour 
obtained from the significance map is shown in the lower 
plot Fig. ITHb ') . The arrows are the estimated angular reso- 
lution (±1(T = ±0.25°). Note that this is larger than that 
of the NGC 253 analysis (0.23°), due to the zenith-angle 
dependence. The Crab was observed at zenith angles of 
around 56°, whereas NGC 253 was observed at around 
6° . The center of the significance map corresponds to the 
Crab pulsar, confirming that our pointing and angular res- 
olution are consistent with our Monte-Carlo simulations. 

Miscellaneous checks were carried out, as described 
here. The agreements between the ON and OFF a dis- 
tributions in a region greater than 30 degrees were com- 
pared with the observation times listed in Table. |21 The 



result is 



dividual observation was calculated 



1.02 . The signal rate for each in- 
and the results are 



plotted in Fig. E| They are consistent with the average 
flux described above. 



4.2. Systematic uncertainty in the Energy 
determination 



There were uncertainties in determining the energy scale 
from the total ADC counts. The ADC conversion factor 
was 92^7^, as previously described. The mirror reflectiv- 
ity also had some uncertainty, in both its value (averaged 
over the whole mirror) and its time dependence. A mea- 
sure of the latter could be made from month-by-month 
shower rates. Also, Mie scattering was not taken into ac- 
count in our Monte-Carlo simulations. Considering all of 
these effects, we estimated the systematic uncertainty in 
the energy determination to be within 15% (bin to bin) 
and 20% (overall). These are also consistent with the Crab 
analysis results described previously. 
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Fig. 19. Excesses versus observation times (after cloud 
cut) for all runs. The blank marks were obtained by an 
OFF-source run, and the filled marks by ON-source runs. 
The dashed line is the average flux. 



4.3. Differential Flux 

The energy spectra of TeV gamma-ray sources gener- 
ally have a power-law nature. Therefore, we used the 
logio{Energy) scale in binning events to determine the 
spectrum for NGC 253, rather than energy, itself. The 
energy for each event was assigned as a function of the 
total ADC counts, where the relation between the energy 
and the total ADC counts was obtained from simulations. 
The excess (gamma-ray) events were observed between 0.5 
and 3 TeV. We divided this \ogio{Energy) range using 
equipartition. The number of binnings is 6. In Fig. 1201 the 
derived spectra for both 2000 and 2001 are plotted, and 
are seen to be consistent with each other. 

The systematic uncertainties were estimated as fol- 
lows. The background light sources, such as stars and arti- 
ficial light may have a significant effect on the estimation 
of the differential flux determinations. These backgrounds 
affect each pixels pulse-height distribution. Small contri- 
butions (Poisson distributed) would be added to the sig- 
nal in each pixel. In order to study the significance of this 
effect, we varied the ADC threshold from 300 (default), 
to 350 and 400 (corresponding to 3.3, 3.8, and 4.3 p.e., 
respectively). The differential flux was obtained for each 
case and plotted in Fig. [2] The total excesses of signals 
varied 1652±149 events (300 count threshold), 1429±154 
events (350), and 1034±160 events (400), respectively. The 
acceptances were (as expected) strongly dependent on the 
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Fig. 20. Differential fluxes for the 2000 and 2001 data. 
The blank squares represent data from 2000, and the filled 
squares are from the 2001 data. The triangles are the up- 
per limits. The dotted line is that of Crab nebula for a 
reference. 



threshold value, but the differential fluxes were stable, as 
shown in Fig. |^ 

As a further check that the excess events are due 
to gamma-rays, the following tests were also made. 
We re-calculated the Likelihood-ratios by first adding 
asymmetry, then removing length, and finally, removing 
width. This was done to check whether only one param- 
eter had an unduly large effect on the final signal. The 
resulting fluxes are shown in Fig. |221 A deviation was ob- 
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Fig. 22. Differential fluxes based on various assumptions 
concerning the Likelihood-ratios. 

served when we added asymmetry to the Likelihood-ratio. 
Asymmetry could only be calculated under the assump- 
tion of a point source at the center of the FOV. As we 
had already concluded that the source was extended, the 
fact that the use of asymmetry reduces the significance 
of the signal is not unexpected, and confirms that it is 
not a useful parameter for diffuse radiation. The other 
Likelihood-ratios were in agreement with each other, sup- 
porting our conclusion that the excess is due to diffuse 
gamma-ray emission. A summary of fiux changes in this 
parameter study is listed in Table 0| The mean energies 



Fig. 21. 

fluxes. 



Energy (TeV) 



ADC threshold dependence of the differential 



Table 4. Summary of the flux changes in the parameter 
study at each energy binnings. 



Energy (TeV) 


0.56 


0.75 


1.07 


1.52 


2.19 3.32 


Flux change (%) 


23.4 


12.0 


8.8 


17.4 


72.8 88.5 



listed in this table were obtained by averaging the gener- 
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ated energies of tlie accepted events in the ADC binnings 
in the above-described Monte-Carlo simulation. 

In order to derive the spectrum of NGC 253 self- 
consistently, we adopted the following method. The above 
spectrum was derived by using acceptances derived from 
simulations in which a spectrum was assumed. If we 

fit the fluxes with a differential power-law spectrum, we 
obtained an index of — 3.7±0.3. We then iteratively used 
this value in simulations to re-derive the spectrum. This 
process rapidly converged at an index of — 3.75±0.27. The 
differential fluxes estimated with E^"^-^ input and that 
with iJ-s 'i's j^j.g shown in Fig. |^ Both showed the same 
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Fig. 23. Differential flux estimated by the various energy 
spectra inputs in the Monte-Carlo simulation. The blank 
squares were obtained by E~'^-^, the black squares by 
E~^''^^, and the black circles by the exponential cutoff 
function described in the text. 



reasonable value. We tried to generate events with this 
spectral input to derive the differential flux again. These 
are shown in Fig. EHl (the black circles). The flux deter- 
mination is very stable over a range of assumptions for 
the Monte-Carlo inputs of the energy spectrum. We also 
carried out an iteration with this function, and confirmed 
the convergence to be good. Finally, we selected this to be 
the Monte-Carlo energy spectrum. 

The systematic errors were estimated by varying the 
Likelihood-ratio cut values, as described previously. They 
are listed in Table. |31 These values are larger than those 



Energy (TeV) 


0.52 


0.68 


0.92 


1.23 


1.73 


2.56 


Sys. error (%) 


29.8 


23.4 


34.8 


23.9 


32.6 


67.9 


dE/E (%) 


36.7 


35.2 


33.7 


32.9 


30.6 


32.3 



Table 5. Systematic errors and energy resolutions of each 
energy binning. 



in Table 0] We, therefore, concluded to use these as the 
systematic errors. Also shown are the energy resolutions 
in each bin, which were obtained from simulations on an 
event-by-event basis. These errors are dominated by the 
core distance uncertainties. From here on, the flux errors 
in the figures are the square root of the quadratic sum of 
the statistical and systematic errors. 

The combined flux is shown in Fig. and Table El 
Around 1 TeV, the flux is about one order lower than in 
the Crab nebula, which is indicated by the dotted line in 
Fig. (Tanimori et al.. ilQQSaf) . 



Log ADC region 



< Energy > 
(TeV) ' 



Flux 

(/cmVsec/TeV) 



3.6-3.8 
3.8-4.0 
4.0-4.2 
4.2-4.4 
4.4-4.6 
4.6-4.8 



0.52 
0.68 
0.92 
1.23 
1.73 
2.56 



(3.51±0.72±1.04)xl0"'' 
(1.09±0.34±0.26)xl0"" 
(6.40±1.69±2.23)xl0-i^ 
(2.05±0.78±0.49)xl0-^^ 
(2.65±2.83±0.87)xl0-^^ 
(2.31±8.07±1.57)xl0-" 



Table 6. Energy binnings and differential fluxes. The first 
errors are statistical and the second ones are systematic. 



best fit spectrum. 

An extrapolation of this power-law spectrum to lower 
energies deviates greatly from the measured fluxes and up- 
per limits (Blom et al.. 119991 Sreekumar et al., IX994|; a 
turn-over below the TeV region clearly exists. Physically 
plausible functions exist with a turn-over include spectra 
oc £;-7e--E/i5™<.x and oc E-^e''^/''. Although the former 
is typically used for the spectra of gamma-rays originat- 
ing from 7r° decay, the value of 7 should be greater than 
2.0 according to the present acceleration theories, in con- 
tradiction with the measurements at lower energies. The 
latter form is typical for an Inverse Compton origin. We 
fitted a spectrum of this form with the EGRET upper lim- 
its. The best fit with this function gave a = 0.28 with a 



The goodness of fit for the various spectra were char- 
acterized by the values. The results for various fittings 
are as follows: 



dF , 

- = (2.85 ±0.71) 



X^/DOF = 2.1/4, 



10" 



dE 
a ■ 



ae^/\E/Eo)- 



1.5, 



(E/lTeY) 



fEjh 



(-3.85±0.46) 



[cm"2g-lrpgy-li 



(2) 



-^TeV" 



6 X 10^5 [cm-^s^^TeV^i], £^0 = 0.0002 TeV, 
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Fig. 24. Combined differential fluxes. The dotted line is 
that of Crab nebula observations. The other lines are the 
fitting results. The dashed line is that of a power-law. The 
solid curve is that with an exponential cutoff. 
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Fig. 25. Significance map obtained by this experiment, 
shown by the thick contours. The thin contours are optical 
image by DSS2. 

Survey). The "Significance" is proportional to the inten- 
sity only when the acceptance and the background level 
are uniform in the full FOV. The detection efhciency is de- 
pendent on the offset angle of the assumed source position 
from the centre of the field of view, as shown in Fig. 



b = 0.25 ± 0.01 [VTcV],x^/DOF = 1.8/5, 



(3) 



where the second formula is based on the Inverse Compton 
scattering formula when we constrain the flux at 3/4 of 
the EGRET upper hmit at i;o=0.0002 TeV. Due to this, 
the two parameters in Eq. (Q, a and Eq, were fixed to 
those values. Here we assumed the power law index of the 
incident electron energy spectrum to be 2.0. A better 
was obtained for this fit compared to the single power- 
law fit. In order to explain our flux and that of EGRET 
simultaneously, Eq. Q is one of reasonable choices. 

4.4. Morphology 

The thick contours in Fig. [53 represent the source mor- 
phology obtained from our observations. These contours 
were obtained from the so-called significance map and are, 
therefore, not exactly speaking a morphology. The signif- 
icance map was made from the distribution of the detec- 
tion significance determined at each point, based on the 
assumption that each point in turn was a point-source 
position. The significance was obtained from the differ- 
ence in the a plot (ON- minus OFF-source histogram) 
divided by the statistical errors. The angular resolution 
of this method was estimated to be 0.23° (Icr is a 68% 
confidence level). The telescope pointed at the center of 
NGC 253. Also shown by the thin contours is an optical 
image obtained by DSS2 (second version of Digital Sky 
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Fig. 26. Efficiency versus offset angle of the gamma-ray 
source position from the center of the field of view. 

We checked the effect of the background light. The op- 
tical magnitude of NGC 253 has magnitudes of ms ~ 8 
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and niyis ~ 7.1. Even if it was concentrated on one point, 
the background level due to this was lower than our sen- 
sitivity (Stars fainter than a magnitude of 5 could not be 
detected in either the scaler or ADC data). We also note 
that the lower cut on Distance was 0.5°, and so pixels 
around the center of NGC 253 were generally not used 
for the analysis. When observing NGC 253, the bright- 
est stars in the FOV have magnitudes of 5.6. However, 
we observed some effects from a group of faint stars (each 
of magnitude ~6) in observations of another target, which 
deformed the shapes of the parameter distributions. These 
effects are believed to be removed by the hot pixel rejec- 
tion algorithm, described in section [3. 61 This was demon- 
strated in an analysis of the Crab nebula data. Although 
the visual magnitude of the Crab nebula is 8.4, a bright 
star (magnitude 3.1) is located within the FOV of our cam- 
era. Despite this, we were able to derive a significance map 
consistent with the other measurements fFig. ll8|) . Because 
we cannot rule out the possibility that "hot pixels" may 
deform the significance map, we can not definitely derive 
the morphology of the gamma-ray emitting regions from 
observations with only a single telescope. 

Fig.|2Z|shows the acceptance of gamma-ray-like events 
as a function of the Distance upper cut values (minimum 
cut is 0.5°). From this figure, we tried to estimate the 
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Fig. 27. Signal yield as a function of the upper range of 
the Distance cut. The yield was normalized to the total ex- 
cess. The points with error bars were obtained by this ex- 
periment. The solid line was obtained from a Monte-Carlo 
simulation with a point source assumption. The dashed 
line is the case for a "diffusion angle" of 0.4°. 



spatial extent of the gamma-ray emission. We proceeded 
by fixing the minimum cut value of Distance at 0.5° (as 
used in previous Whipple and CANGAROO analyses). 
The maximum cut value was varied between 0.6 and 1.5, 
and the excess events for each cut were plotted by the 
points with error bars. 

We checked the consistency between the experimen- 
tal data and the source diffusion assumptions. At first, 
the solid line in Fig. |23 was obtained by a Monte-Carlo 
simulation with a point-source assumption. The observed 
distribution is clearly broader than this. 

From the significance map (Fig. I25|l . the correlations 
between the orientations of the TeV emission and optical 
image were calculated and the standard deviation of the 
long axis was 0.37° and of the short axis was 0.24°, re- 
spectively. The long axis was inclined by -1-30° from the 
horizontal axis. This is slightly larger than that of the op- 
tical image; however, we do not believe that the difference 
is significant. An analysis of the map of the number ex- 
cess events yielded similar results. We then carried out 
a Monte-Carlo simulation based on various assumptions. 
We varied the extent of the emitting region by smearing 
the gamma-rays' incident angle with a Gaussian. Our data 
are consistent with "diffusion angles" of between 0.3-0.6°. 
The case of 0.4°, shown in Fig. EH by the dashed line, is 
consistent with the observation. However, it is necessary 
to wait for future stereo observations (CANGAROO-III 
and H.E.S.S.) before any quantitative estimates of the ex- 
tent of the emission region can be made and the mor- 
phology studied. We can concluded that the emission is 
consistent with Gaussians between 0.3-0.6 °, which cor- 
respond to 13-26 kpc at a distance of 2.5 Mpc. Changing 
the input gamma-ray's spatial distribution from Gaussian 
to a rectangular shape gave similar results, i.e., Fig. |57| 
can be reproduced by a 0.4° rectangular spatial emission. 
This can be understood by the efficiency reduction due 
to the limited FOV. The acceptance for our telescope was 
reduced from an offset angle of 0.5°. The electrons of GeV 
energy associated with NGC 253 was reported by radio 
observations of Hummel, Smith, and van der Hulst II1984II 



and Carilli et al. I|1992(l . The size of the emission region is 
similar to our result. An interpretation of this results can 
be found in Itoh et al. (|2003,). 

5. Conclusions 

In this paper we have concentrated on technical details 
concerning our observations and analysis. Statistically sig- 
nificant signals of gamma-rays from the nearby starburst 
galaxy NGC 253 have been detected. The differential flux 
shows a turnover below 0.5 TeV. The spatial distribu- 
tion of the gamma-ray emission is broader than that of 
a point-source. This is consistent with a width of 0.3-0.6°, 
corresponding to 13-26 kpc at the location of NGC 253. 
A more detailed physical interpretation is presented else- 
where (Itoh et al., 12003 ■ 
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